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Abstract 

The loop structure of two point Green's functions is investigated in the Wess- 
^vq \ Zumino model in the formalism where the auxiliary fields are integrated out. In 

the usual frame of perturbation theory the deviation from the non-renormalization 
theorem is explicitly shown. It is shown that Ward identity are not satisfied in 
this approach. Further we go beyond perturbation theory by solving a system 
of regularized Schwinger-Dyson equations (SDEs). The mass splitting between 
fermions and bosons, which was already observed in perturbation theory level, is 
' further enhanced. 
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^ : 1 Introduction 

o - 

In this paper, we investigate the loop structure of the field theory model in- 
troduced by Wess and Zumino [lj in formalism where the auxiliary fields are 
integrated out. Our aim is to explicitly check the (m-) equivalence of quantum 
theories defined with and without the use of auxiliary fields. The later approach 
(denoted as II in the next text) is based on the quantization of the Lagrangian 
model which has been obtained by eliminating auxiliary fields F ) G within the 
help of their classical equation of motion. The first is the well known manifest 
supersymmetry (Susy) invariant (denoted as I in the next text), which in its own 
has been examined many times for different reasons. 

First let us ignore quantum correction and summarize the classical theory. In 
such case, the classical action is in our notation [2] 



&c = -${dA) 2 - $(dB) 2 - + + 
+ m(FA + GB - i^>) 

+ g(FA 2 - FB 2 + 2GAB - ^A - i^j^B) , (1) 

where classical Lagrangian involves a scalar field A, a pseudoscalar field 5, an 
anti-commuting Majorana spinor ^, and two auxiliary scalar F and pseudoscalar 
field G both of canonical dimension mass 2 . Under supersymmetry transformation 
each line in Lagrangian jC^ transforms as a total derivative, composing thus Sq 
from three independent Susy invariant terms. 
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Using a classical field equation of motion for F, G one can obtain the following 
classical action 

S£ 7 = / d A xCf II 



^ " "2 



1 



d^Ad^A + d^Bd^B + jfo/j + m 2 A 2 + m 2 B 2 + 
mgA(A 2 + £ 2 ) - ^-(A 2 + £ 2 ) 2 - - ig^j^B . (2) 



It is widely believed that theories I and II should describe the same physics even 
when they are quantized, i.e. it is assumed that the equations of motion for the 
fields F, G survive quantization. 

In the paper [H| the classical theory I was quantized via path integral and 
considered generating functional for the Green's functions is 

Z^J 1 } = J VAVBVFVGV^jexp^ s Q , (3) 

where the quantum action is defined as usually 

= S T C + J d 4 x[AJ A + BJ B + FJ F + GJ G + #] , (4) 

where we have introduced space time dependent function J 7 (x) = {Ja^Jb^Jf^Jg^v) 
associated with the field $ = (A, B ) F. G. ifj), where J 1 includes c-number exter- 
nal sources of bosons and rj is an anticommuting Majorana spinor. Imposing the 
transformation rules for J 1 then the Susy can be retained and the well known 
Ward identities can be write down (for Susy Ward identities see original paper 

m 

Generating functional for the model II is given by 



Z n [J n ] = (0 



Texpli J d 4 x [J A .A + J B .B - fjip] J 



0) , (5) 



where in contrary to the previous case , the Lagrangian of the model II has 
the ordinary degrees of freedom, i.e the all fields have the usual kinetic terms. 
As consequence of elimination F ) G the all fields have now non-zero canonical 
conjugated variable with non-vanishing equal-time canonical (anti) commutation 
relations. From this, using a standard treatment one can arrive to the alternative 
expression for the functional © 

Z n [J n ] = ^JvADBTtyexpliSg} 

S% = S'J + Jd 4 x[J A .A + J B .B-fjiP] (6) 

where Jn = (Ja, Jb,v) are classical sources of Wess-Zumino physical multiplet 

<j>=(A,B,il>). 



1 Introduction 



3 



Assuming no anomaly in I, the Ward identities 

n A (p 2 ) = n B (p 2 )=p 2 r 1 ( P 2 ) (7) 

should be respected by the renormalization. In general these identities constrain 
various renormalization constants. Furthermore many of these constants are not 
actually necessary. Already in the paper [3J the non-renormalization theorem 
(NT) for WZM was proposed. This NT simply state that renormalized and bare 
quantities are related as 

3 

rn r = Zrri] g r — Z 2 g (8) 

where Z is the field strength renormalization constant which is common to all 
field, i.e. i/; = Z l / 2 ip r and the same is valid for the fields A and B; m r is 
renormalized mass and g r is the renormalized coupling constant. In the words 
the equations ((El) means that all diagrams belonging to vertices should be finite 
and the propagators could be divergent in very uniform way. In this place we 
should stress that the NT, ie. the actual finiteness of vertex functions, are proofed 
by an explicit evaluation of diagrams in superspace ^fl (see also [H] for 
an indication and for instance [Zj for the renormalization without supergraph 
technique) i.e. the appropriate proof is unavoidable based on the auxiliary field 
formalism represented by I. 

Up to my knowledge a proper examination of infinities remains undertaken 
in Susy field theories without the explicit help of auxiliary field formulation. 
However, dealing with theory with true physical degrees of freedom pertained, 
i.e. with the functional II, we are not able to reproduce the NT at all. In the next 
section we will demonstrate the deviation from the WTI and the NT explicitly 
already at one loop. 

Although, we observe that one common Z is still sufficient to renormalize all 
field content of WZM, the explicit calculation suggests that the masses of boson 
and fermion sector necessarily differ due to the quantum corrections. Explicitly, 
at one loop level we need to introduce the boson mass renormalization function 
Z mb , such that 

m b = ^m b r ,; b = A,B (9) 

where = mf are renormalized boson masses when Z mA = Z mB is naturally 
taken. Remind here for completeness, that the relation = m^/Z is still 
satisfied since the field F, G do not couple to In order to renormalize the 
WZM II, we must add the counterterms which break the original supersymmetry 
invariance. Susy, if formulated by the quantum action in II, does not survive 
quantization . 

The paper is organized as follows. In the next section we write down the one 
loop perturbative result for selfenergies. In the section 3 we go beyond the frame 
of perturbation theory and derive the set of Schwinger-Dyson equations (SDEs). 
After a simple truncation of SDEs we use some regularization techniques to solve 
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the SDEs numerically. Proposed method offers a simple estimate of Green's 
functions behavior without loosing connection with the 'bare' Lagrangian. 



2 1-loop corrections 

For the discussion of loop structure and the renormalization it is useful to define 
generating functional of one particle irreducible (1-PIR) Green's functions. This 
is obtained from the generating functional W of connected Green function via 
Legendre transformation 



T[R] = W[J]- J d 4 x(J A .A + J B .B 

W[J] = -ilogZ[J], (10) 

where R are the semiclassical vacuum expectation of cp in the presence of J already 
defined in the introduction. They satisfy 

_ sw sw sw 

- jt a ' Rb -jt b ^-~^> 

j - — j - v - (12) 
S Ra 8 Rb S R^ c 

First of all, let us mention that in the formulation I fermions does not interact 
with auxiliary field. Therefore there is no difference when compare one loop 
fermion selfenergy in II and in I (starting from (fTUj) we will omitted index II 
since from now we will deal only with the theory II and no confusion threats ). 

In addition we start with the calculation of the propagator Gb- The calcula- 
tion of Ga will proceed similar steps. (We will relegate the appropriate derivation 
of the propagator Ga in fully nonperturbative fashion to the next Section. The 
perturbation result then easily follow). Before doing this explicitly we conven- 
tionally define the proper two point function 



Y B {x-y) = 



S 2 T[R] 



5R B (x)5R B (y) 



j=o 



d 4 p 



e ip-(*-y)Y B ( p 



2\ 



T B (p 2 ) = T^ s -U B (p)=T^ s + Tf op + ... 

T das = p 2_ m 2^ ( 13 ) 

(and the same we define the propagator of A) noting that the identities G B 1 (p) = 
^ b{p) )G~a {p) = Ta(p) are exactly valid even in the presence of external sources. 
The fermion propagator is defined as 

S- 1 (p) = f (p) + Tf (p) = - m - E(p) , (14) 
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where I\, T 2 are two independent function and E(p) is the conventional selfenergy. 
Equivalently, it is sometimes convenient to introduce the following formula 

Hp) 



Sip) 



M(p) 



(15) 



where the so called renormalization wave function is defined through the eq.: 
T~ x = r« and mass function is M = T^/T^. 

The appropriate contribution to Ub can be represented as 



+ 



+ 



+ 



» — ' 



o 



-/ 



d A i 

W)' 



Gg 2 



I 2 — m 2 + is 
Am 2 g 2 



+ 



2g 2 



I 2 — m 2 + is 



{{I 



p) 2 



is) (I 2 



+ 



■ is) 
Ag 2 p 2 



/ _ ro2+K((i _ p)2 _ ro2+fe)(P _ m2+I£)JJ . («) 

where each term in the each bracket [] in $TBJ) is associated with given Feynman 
diagram depicted in the first line (in given order). The dashed (dot) line corre- 
sponds with the propagator of A(B) field, the solid line stands for the fermion 
propagator. Convenient shift and some obvious algebra has been done ( we clear 
out the numerator l.p by the use of the identity 21. p = G~ l {l) — G~ l {l — p) +p 2 ) 
in expression for the fermion loop. The appropriate cancellation of quadratic 
divergences between scalar tadpoles and fermion loop is obvious. Summing all 
together, the complete one loop result is given by the following logarithmically 
divergent expression 



n B (p) 
Hp 2 ) 



a 



-0 2 + m 2 )l(p) 

7T 

rdH 
-i / - 

J 7T 



(17) 



1 



2 {{I 

A 2 



p) 2 — m 2 + is) (I 2 — m 2 + is) 



Am 2 



lim / duu- 

a^oo j4 m 2 — cj + ie 



(18) 



where we have defined a = g 2 /(47r). 

For completeness we also review the appropriate set of relevant Feynman 
diagrams which complete one loop contribution to 11^. These are 



n A (p 2 )= — — + — -- — + 



--o 
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• * • 

+ ;- + ... 

^ • . • 



(19) 



The evaluation of Ha simply gives 



U A (p 2 ) = Il B (p 2 ) . 



(20) 



And recall that in our notation the fermion selfenergy can be evaluated as 



In the limit A — > oc in ([T%|) the infinite term m 2 ^I(p 2 ) in the expression 
ffTHjl for boson selfenergy can be removed only by renormalization of the masses 
m A and m B in the WZM. The NT as well as Susy WI are not valid for theory 
II. Assuming g is a rather small number, i.e. g << 1 then the perturbation 
theory inspection provides the strong evidence (if not proof) of non-equivalent 
formulation of the supersymmetric quantum field theory defined throughout I 
and II. Note, it does not mean that one of the the theories is incorrect. The 
result should be read as that the quantum field theory II is not equivalent to 
the theory I wherein the classical equation of motion for phantom fields F ) G are 
simultaneously satisfied (as it was supposed in [T]). 

3 Masses beyond perturbation theory 

In this Section we derive Schwinger-Dyson equations (SDEs) for the propagators 
of the WZM II. These are the first of infinite system of quantum equations 
of motion which relate Green's functions of theory with themselves. If solved 
exactly, we would obtain the full information about the theory. The technical 
impossibility of such a procedure is obvious and we truncate the system of SDEs 
in a way that if we approximate the exact Green's functions inside the loops 
by the bare ones then the one-loop results obtained in the previous Section are 
reproduced. 

Before presenting the derivation we anticipate the resulting equations graph- 
ically. The set of SDEs for two point Green's functions can be diagrammatically 
represented as: 



np)=i>-i{p 2 )- 



(21) 



G- A 1 (p 2 )=p 2 -m 2 -U A (p 2 ), 




3 Masses beyond perturbation theory 



7 




+ —( Q 0-- | — { Q 0.- 



(22) 



n B (p 2 ) 



l 



+ . . . . .v. 



/•••. \ 



+ ...a. — + ...( 0. 



(23) 



5" 1 (p)=^-m-E(p) 



S(p) 



(24) 



where now (in contrary to the previous Section) the all internal lines inside the 
loops represent the fully dressed propagators. The diamonds and bloops stand 
for dressed quartic and trilinear vertices which satisfy their own SDE. The above 
mentioned truncation of SDEs simply means that in the equations . . 
the two-loop diagrams should be neglected and dressed vertices should be replaced 
by the classical ones. 

The derivation of SDEs is tedious but necessary task we now turn. For this 
purpose we accommodate the elegant and waterproof functional method originally 
developed by Symanzik and others [H] . Review of this formalism can be found 
in some standard introductory textbook (e.g. jHj). The SDEs can be derived by 
further variations of the following compact formula: 

ST[R] _ SS C 
5R^(x) 5(f>i(x) 

where is operator acting on the right 

S-(x) = R>(x)-i U\ d ' W[J l — 



(25) 
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where J^. is an associated external source of R^. is the classical expectation 
value of In words, the right handed differential operator <^ replaces the 
Heisenberg field operators (f) positioned on the left in the classical equation of 
motion. In order to obtain the physical Green's functions the sources are switched 
off at the end of the calculation. 

Now we are going to explain some details for the equation for G^ 1 . For the 
case of space-like derivatives of (j) = A we can use per-partes integration as usually 



ST[R] 



+ m 2 )R A (x) + g^(x)^(x) (27) 



SR A (x) 

+ mg(3A 2 (x) + B 2 (x)) + 2g 2 A(x)(A 2 (x) + B 2 (x 

A( X ) = R^)-ij^ sj ™lf A(y) J A(v y 

A, , „ , , f ,1 PW[J] S 

B( X ) = *B(*)-iJdy My)SRB(v) , 

where we have used obvious identity 

${x)l = R 4 ,{x). (28) 

In addition we will use convenient shorthand notation for the measures and 
dirac delta in Minkowski spacetime. 

dx = d 4 x; dp=— — ; 5 xy = 5 4 (x - y). (29) 
(2tt) 4 

Also we introduce the abbreviation for the variations of the generating function- 
al: 

s n r[R] = /-',,. 

6R il (x 1 )...6R in (x n ) ' 

snw[J] - (30) 



5J 4>ii {x 1 )...8J 4>in (x n ) 
Once more variation of the Eq. pTj) with respect to Ra gives 



r? y f) A = -d^S xy - m 2 5 xy - U A (x, y) (31) 

s 

Ka(x, y) = ^-y {mg[3A 2 (x) + B 2 (x)] 

+2g 2 A(x)[A 2 (x) + B 2 {xj\ + g$(x)4>(x)} (32) 
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Switching the external sources J the equation (|3T|) represents the Fourier 
transformation of the inverse propagator G^ 1 

G?(x, y) = rf^| fl0=o = / dpe*< x -y) [p 2 - m 2 - U A (p)} , (33) 
since in our notation 

K% U=o = " / dpG A ( P )e ip ^ . (34) 

Recall at this place general functional identity 

/ dzFf^W^ = -S xy (35) 

which is valid in the case of the presence of the external sources J to. 

The cubic A 3 term in the Lagrangian is responsible for the following contri- 
bution to Ha (i.e. the first term in 



8A 2 {x) 
SR A (y) 



l 



SR A (y) 

i 



R 2 A (x)-iW ( A x 



(x,x) 



= 2R A (x)S xy 
The second required variation in (pTTj) is 



(36) 



8B 2 (x) 
SR A (y) 



l 



SR A (y) 



R 2 B (x) 



dzduW^ z) Tf y fX B W^ 



(37) 



where we have repeatedly used the relation 



SWf 1 * 2 v 

{X1,X2) 

sn( y ) 



(38) 



which is valid for any Q, however thorough this paper Q always stands for some 
R. Note that eq. (p$£|) simply follows from the identity (|3Kjl . Taking J, f? = 
then the obtained integrals in eqs. (|37>|) .(j3*T j) correspond with the third and the 
fourth diagram in the diagrammatic expression for selfenergy 
The terms which follow from the quartic interactions are 



SA 3 (x) 
SR A (y) 

+H) 2 



1 



SR A (y) 



R A (x) - ZiR A (x)Wfe 



AA 

x) 



(39) 



MzduW*%W*%Tf*X A K% 
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+3H) 2 j drdsMzduW ( i A r) r«^ 



SA(x)B 2 (x) l 
SR A (y) 



SR A (y) 



5 



R A (x)R%(x) - iR A (x)W; 



(40) 



+H) 2 / dtdzduW { i%w^%rf t -X B K% 

= R 2 B (x)S xy - iW* B x) 5 xy 

-iR B (x) J dzduW? x %Tf y f^W™ x) 

+H) 2 / MzduW^W^Y*^ 

+2H) 2 J drdsdfdzd^ 

+ H) 2 / drdsMzduW { i%r*f^ . 



Setting R to zero then the first line after the second equality of the above 
expressions (|39j) . (JlOj) is responsible for generation of the tadpole diagrams wherein 
the vertex is undressed. The third lines of the above expressions correspond 
with the sunset diagrams wherein one quartic vertex is dressed. Considering the 
prefactors the sunset diagrams contribution to n(x,y) is as following 



where T^^a 2 b 2 are the exact proper quartic vertices. The momentum cor- 
responding to external leg of the full vertex is ingoing, while all the momenta 
associated with lines inside the loop are outgoing from this vertex. 

The three remaining contribution in corresponds with twoloop di- 

agram of the same topology. The last term in (JSHl) corresponds with the last 
diagram in (22), while the both remaining terms in (^OJ) correspond to one to last 
diagram in rel. (j22j). 

The last term in eq. corresponds with the fermion loop diagram depicted 
in ((221) • Suppressing dirac indices the appropriate contribution to Ha reads 



(i) 2 2g 2 J dtdqdpj p - {x - v) G A (p ~ l)G A {l - q)G A (q)T A i (p; q ,l- q ,p-l) 
{i) 2 2g 2 J dldqdpe ip -^G B {p - l)G B {l - q)G A {q)Y A , B ,{p- q,l- q ,p-l) 



iTr JdzduT^Wt^Wl^. 



(41) 



5A(y) 
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The derivation of higher vertices proceed through the further differentiation 
with respect of a given i?'s. However in this paper we need to know only clas- 
sical (tree-level) vertices, zero order of the Planck constant expansion result 
must correspond with the variation of the classical action with respect to the 

a given original fields, i.e. A,B,i/j. Using the obvious and convenient nota- 

Ra Rb 

tion r A k B n (xi...x n +k~) — rV Bn '\ we list the appropriate vertices here 
for completeness: 

T A s(x, y, z) = -6mg5 xy S yz ] T AB 2(x, y, z) = -2mg8 xy 8 yz 

T^ A (x, y, z) = -2gS xy S yz ; T^ B (x, y, z) = -2gS xy S yz ; 

r A *(x,y,z,t) = -12g 2 5 xy 5y Z 6 zt ] T A 2 B 2(x ) y, z, t) = -4:g 2 5 xy Sy Z 6 zt ] 

T B 4(x,y,z,t) = -12g 2 5 xy 5 yz 5 zt . (42) 

Substituting the vertices (^2j) into their right places, we can write down the result- 
ing expression for the function II^. For practical purposes we Fourier transform 
the result into the momentum space. After the little algebra we get 

U 4 P) = -ri- fdlU,f (l) f}± - ^1)1-* + (M(l) + MW\ 



l 2 -M 2 {l) (I 2 - M 2 (l))(q 2 - M 2 (q)) 

+6G A (l) + 2G B (l) + lSm 2 G A (l)G A (l -p) + 2m 2 G B (l)G B (l - p)} , (43) 

where the second line follows from bosonic loops while the first line follows from 
the fermionic loop and where we used notation q = p — I. 

The remaining SDEs can be derived from (|25J) by very similar fashion as the 
equation for G A . From the following definition 

5 2 Y[R] 



dR B (y)dR B {x) R=Q 

we can obtain the equation for the selfenergy function H B . In our approximation 
of SDEs the result reads 

G B \p) = p 2 -m 2 -n B (p) 

n B (p) = -ig 2 J dl{QG B (l) + 2G A (l) + Am 2 G A (l)G B {l - p) 

4 ' w v /L . r„,:r: vv , ;H 2 ^ } (45) 



, mm _ , hi)h<i)[-p 2 + mi) - M(g)r 

'l 2 -M 2 {V) (I 2 - M 2 (l))(q 2 - M 2 (q)) 



The majorana fermion mass function Ai as well as the renormalization wave- 
function T in the bare vertex approximation of the fermion SDE are 

noting that there are only one loop diagrams even in the exact case. 
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4 Numerical results for SDEs of WZM 



In this section we present numerical results for the propagators as they have been 
obtained from the solution of the SDEs. To obtain this, we convert SDEs from 
Minkowski space into the Euclidean space via performing Wick rotation. The 
resulting set of rotated equations is presented in the Appendix A for a reader's 
convenience. Here, we should mention about possibility to solve SDEs directly in 
Minkowski space [TO], however this far nontrivial approach is until now developed 
for renormalized theory only. Therefore at this stage of our calculation we deal 
with Euclidean formalism, which method leaves us with the solutions for space- 
like momenta only. 

Before describing a quantitative face of the model we anticipate basic re- 
sults. The most striking fact is the unexpected behavior of the boson selfenergy 
functions Ha,b- In accordance of our expectation they are degenerated, i.e. the 
equation 11^ = Ub is valid (with negligibly small numerical error). On the other 
the boson selfenergy appears to be fairly unrelated with the fermion one. The 
already perturbative observed violation from Susy Ward identity is further en- 
hanced. Again we will present the solution for unrenormalized but regularized 
Green functions. (Notation: solving the renormalized set of SDEs one can always 
approximately preserve Ward identity at least in the vicinity of the renormaliza- 
tion. This is not case we are interested in.) 

Susy invariant regularization of I was already proposed in the paper The 
regulator Lagrangian term of the form 



has been added into the original Lagrangian I ( Ckin corresponds with the first 
line of QJ. This approach -albeit useful in perturbation theory analyzes in the 
limit A — > oc- becomes extremely inconvenient in nonperturbative study like here, 
especially when A remains finite. In fact, due to the number of time differentia- 
tion, new degrees of freedom appear and it is rather nontrivial task to recognize 
particle content of such theory ( associated with the S-matrix singularity rather 
then the number of the fields in the Lagrangian). However, the theory II does not 
respect Susy due to its radiative corrections breaking. We are not so restricted 
by the choice of the regularization scheme and we use two very simple regulator 
functions in our set of Euclidean SDEs. As a first we have explored hard cutoff, 
i.e. the momentum integrals in foUl become 



As a second prescription we use the smooth regulator function, such that 



A 2 



(47) 




(48) 




(49) 
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Figure 1: Test of numerical precision. Square root of the boson selfenergy 
^a,b(p 2 ) as has been obtained with various number of mesh-points, i.e. (num- 
ber of discretized momentum points ^number of discretized point of cos 9) Here 
A 2 = 10 6 ra 2 . 



\XjZ are the squares of Euclidean momenta defined below eqs. (|5()jl ). 

Up to a few percentages accuracy, we have found, these two regularization 
schemes become identical for A >> m. The presented numerical results were 
calculated within the use of second approach, i.e. with regulator (flHjl . 

To integrate the equations numerically we use Gaussian quadrature method. 
In the case of usage of regulator function the largest momentum was typically 
p >> A (while the obtained functions have physical only in the range p < A). The 
equations (|5Uj) have been solved by the method of iterations. In this treatment 
we need a knowledge about the functions on a grid which differ from the one 
we choose to discretize of momenta y — > yi. For this purpose we interpolate the 
appropriate functions inside the kernels of SDEs. The numerically convergence 
was observed already for relatively small number of mesh points. This is shown 
in the figure The renormalization wave function T is displayed in the figure 
noting that due to the degeneracy of scalar and pseudoscalar the mass function 
is fully given by its inverse i.e.A^ = mj T . The boson selfenergy is shown in the 
figure For a very large cutoff it behaves like a constant in the low energy 
where s = p 2 « A 2 . The results were calculated for rather small coupling 
a = 0.01 and for three different values of cutoff A 2 /m 2 = 10 5 , 10 10 , 10 20 . 
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Figure 2: Renormalization wave function T for two values of the coupling con- 
stant a = and ultraviolet cutoff A 2 = 10 8 m 2 . m is the Lagrangian mass. 




le+40 



s 



Figure 3: Square root of the boson selfenergy YIa,b(p 2 ) for & = 0.01. Figure 
shows the dependence on the ultraviolet cutoff A. The result for a = 0.1 and 
A 2 = 10 8 m 2 is added or comparison. The results have physical meaning only 
bellow the cutoff ,i.e. on the left side from the maximum of a given line. 
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5 Conclusion and speculations 

We studied quantum corrections to the propagators of Susy model without the 
explicit presence of auxiliary fields. In the section 3 we have shown that the 
next to leading corrections correction lead to different prediction then known 
from manifest Susy invariant treatment theory I. Calculated Green's functions 
do not respect Susy Ward identity own to the formalism I. The origin of this 
can be easily traced: the preferable choice of unperturbed Lagrangian [I] does not 
correspond with with free Lagrangian II. 

In addition we have shown that this 'unexpected discrepancy' between I and 
II (however it was somehow indirectly intimate already in the paper [I]) be- 
comes even more crucial when the solution of SDEs is considered. Of course, the 
appropriate solution is only certain approximation of the full nonperturbative 
solution (e.g. it is truncation of SDEs system dependent), however, in agreement 
with experience from another quantum field models, there is no large chance for 
some significant changes when further improvement will be made. Some lattice 
simulation, recently undone, but hopefully made in the near future could out- 
shine shades at this area (for some progress see the paper [II]). The observed 
violation of non-renormalization theorem has further non-trivial consequences. 
This is because the coupling constants of cubic interactions (in II) are propor- 
tional to bare Lagrangian mas m. These need to be specified when dealing 
with renormalization. The already mentioned lattice formulation of the problem 
could be actually credible guide. At this place we should mention the paper 
|12| where the authors observed the splitting of the renormalized coupling con- 
stant as g — > Qaaf ? 9ag b 2 ? 9Aipip (remind, the authors of this paper dealt with I). 
This is an another indication that a naive formal quantization and associated 
functional method applied on auxiliary fields (i.e. the field without canonical 
conjugate momenta) can fail in prediction when compared to the result obtained 
from the sophisticated first principle method (here directly from the discretized 
path integral). 

Although the model studied here has obviously no phenomenological interest, 
it has a large historical impact on a part of particle physics. At given time an 
extreme effort and many speculations have been pursued to uncover and explain 
possible mechanism of Susy breaking O E3J ESI EE E| ?• • • • It seems very likely 
that such troubles can be overcome if one use honestly constrain ,i.e. equation 
of motions for auxiliary field content, before the quantization and thus then deal 
with the usual quantum field theory expanded around the classical vacuum, where 
the nearest excitation are just free physical particles. 
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6 Appendix A 

SDEs for WZM in Euclidean space 

Converting (|43|45l46j) to the Euclidean space (7 = iUe) we should obtain 

q, poo pw q> poo pit 

H-A(y) = n — n / dx d6xsm 2 6J A ; U B (y) = — — / dx d6xsm 2 9J B 

Z7T JO JO l-K z Jo JO 

J A = ^QG A (x) + 2G B (x) - 18m 2 G A (x)G A (z) - 2m 2 G B (x)G B (z) 



F(x)F(z) A T{x)T{z) [y + (M (x) + M (z) f 



+ 4 



x + M 2 (x) (x + M 2 {x))(z + M 2 {z)) 



j (50) 



Jb = [qG b (x) + 2G A (x) - 4m 2 G A (x)G B (z) 

o W) , Hx)Hz)[y + (M(x)-M(z)) 2 } \ 

x + M 2 (x) (x + M 2 (x))(z + M 2 (z)) J 1 ' 



where z = y — 2yJxycos6 + x 5 = [x + m 2 + Ha,b{%)}~ 1 
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